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Abstract 

We develop a generalized version of heavy-baryon chiral perturbation theory to describe pion-nucleon 
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scattering in a kinematic domain that extends continuously from threshold to the delta-isobar peak. The 
P-wave phase shifts are used to illustrate this framework. We also compare our approach with those in the 



literature that concern the delta resonance. 
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I. INTRODUCTION 



Pion exchange provides the long-range components of nuclear forces, and crucial to its under- 
standing is pion-nucleon (ttN) scattering. A prominent feature of the latter is the delta resonance, 
A(1232), a peak in the elastic cross section at the center-of-mass (CM) energy ttia = rriN+5 ~ 1230 
MeV, where 5 ~ 290 MeV is the nucleon-delta mass splitting [1]. Our goal here is to investigate 
ttN scattering from threshold up to the delta resonance in an effective field theory (EFT). 

It is well known that resonances can be studied by considering the unitarity and analyticity of 
the S matrix. Assuming for simplicity that there is only one decay channel for a resonance, when 
the CM energy E is close enough to the resonance, the T-matrix element in the resonance channel 
can be written in the form of a Breit-Wigner formula plus a non-resonant background (see, e.g., 
Refs. flfl), 

T(E) = _i {exp [2i e(E)) - 1} = - E _ E T R + iT/2 (1 + iTs) + T B , (1) 



where 9(E) is the phase shift, Er and T/2 are real numbers that represent the energy and half- 
width of the resonance (that is, the pole position in the complex-energy plane), and Tr is a 
complex number often called the non-resonant background amplitude, which can be written as 
Tb = — i[exp(2i#o) — 1] i n terms of a non-resonant phase shift 9q. 

Though based on general principles, the direct application of Eq. ([T]) relies on a few assump- 
tions. First, in many cases where En and T are unknown a priori, they are free parameters of 
Eq. (pQ) despite the fact that they are related in some underlying theory (which is Quantum Chro- 
modynamics (QCD) in the delta case). In the absence of Tb, one could precisely determine Er by 
where the phase shift passes through ir/2 (hence a sharp peak in the cross section). However, Tb is 
in general non- vanishing so one has to make an educated guess about where the resonance window 
is centered. Therefore, the values of Er, T, and Tb based on Eq. ([T]) may be model-dependent. 
Second, one has to assume how close to the resonance is "close enough". Particularly, is Er small 
enough so that Eq. ([TJ is valid even at threshold? If Eq. ([1]) works only in the resonance region, 
one has to assume that the data points one wants to fit are sufficiently close to the resonance. Or 
can we extend Eq. ([1]) simply by allowing Tb and/or T to be functions of energy (as it happens in 
general in field theory) and, if so, with what constraints on Tb(E) and T(E)7 If T depends on the 
energy, then the position of the pole is not given by ER — iT(ER)/2, and the physical interpretations 
of Er and T(Er) are unclear. In the case of the delta, one can show explicitly [4] that Er and 
T(Er) are indeed unphysical in the sense that they depend at two- loop level on the choice of fields. 
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Given that one cannot, at present, straightforwardly solve QCD at low energies due to the 
difficulties posed by the large strong-coupling constant and the small pion mass m^, EFT is a good 
alternative for describing low-energy nuclear and hadronic physics consistently. Following several 
seminal papers EFTs have been developed as model-independent approximations to low-energy 
strong interactions, which can be systematically improved by a series in powers of Q/Mqcdi where 
Q refers generically to small external momenta and Mqcd ~ 1 GeV is the characteristic QCD 
scale. The loyalty of EFTs to QCD, the underlying theory of low-energy strong interactions, is 
manifested by the fact that EFTs inherit all the symmetries of QCD, among which chiral symmetry 
is probably the most nontrivial. For reviews, see, for example, Refs. 0, Qj. 

A particular EFT, chiral perturbation theory (ChPT), specializes, and has proven quite suc- 
cessful, in processes involving at most one nucleon jf|. At energies E close to the irN threshold, 
E ~ rriTT, that is, for Q around the pion mass, the delta dynamics can be considered short-range 



wit 



i only pion and 



physics, which amounts to treating 5 as a large scale, Q ~ m n >C 5. ChP r 
nucleon fields has been extensively applied to irN scattering near threshold 8Hlll|. However, it is 
a perturbative expansion in powers of Q/5 and which should converge slowly because in the 

real world 5 ~ 2m n . The slow convergence then contaminates pion-exchange nuclear forces [l^ . 
where pion energies uj are small compared to the irN threshold energy, oj <C m^. 

One can improve convergence by considering the delta as an explicit degree of freedom, in 
which case one can take Q ~ ~ 5. One can show that hadronic scattering amplitudes can be 
written as perturbative series in powers of Q/Mqcb, "Itt/^qcd, and 5 /Mq cr> Q, Q], and that 
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ldj . The positive role of 



pion-exchange nuclear forces display a good convergence pattern 

an explicit delta field in irN scattering within ChPT has been explored fl?) ] and demonstrated in 
a fully consistent calculation [181 ]. 

Nevertheless, this perturbative expansion diverges around the resonance, where the delta goes 
on-shell. This is not surprising since the perturbative nature of standard ChPT makes it impos- 
sible to describe such a non-perturbative phenomenon. In order to fully describe low-energy irN 
scattering one needs to resum certain terms in the expansion so as to have finite results throughout 
the low-energy region. A non-perturbative treatment of the delta within ChPT was considered in 
Ref. [191 ], where the leading delta self-energy was resummed. However, a systematic resummation 
did not exist until the seminal work of Ref. j^ . where it was justified by a power counting based 



on three separate scales « { « Mi 



netic reactions in the delta region 
221. 



This idea has since been applied to various electromag- 
20 22], but for ttN scattering few results have been published 
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Note that other approaches exist to incorporate the delta (and other resonances) consistently 



with chira 
proach [151 . 



symmetry and field-definition independence. They are reminiscent of the original ap- 



23| to nucleon-nucleon interactions using chiral Lagrangians: a pion-nucleon "kernel" 



is first derived in ChPT to a certain order and then unitarized, for example using the N/D method 
24l ] or the Bethe-Salpeter equation 25]. Power counting is not manifest at the amplitude level, 
but good results for pion-nucleon phase shifts are obtained past the delta region. 

In this paper we realize the delta as a heavy baryon that fulfills Lorentz invariance order by order 
in powers of Q/itin, using the nonrelativistic delta field as one of the building blocks. We employ a 
power counting developed for generic narrow resonances [26], previously ap plie d to the shallowest 

3- I n this power 



P-wave resonance in nucleon-alpha particle scattering at very low energies 
counting we consider only two scales M\ Q ~ 5 ~ and My ~ Mqcd> and expand in Q/M^i and 
Mi /Mhi. Certain contributions are enhanced near the resonance over the standard ChPT counting, 
which continues to apply near threshold. We can calculate the ttN amplitude spanning the two 
regions in a controlled expansion. We illustrate the method by explicitly calculating the first three 
orders of the expansion and comparing with known P-wave phase shifts 28| and scattering volumes 



[29 ] . In the delta region our approach is similar, but not identical, to that in Ref. [201 ] . We compare 
our approach with those in Refs. fill I20I as the differences arise. 

The rest of the paper is organized as follows. In Sec. [II] we describe the heavy-baryon chiral 
Lagrangian that has an explicit delta field. (Some of the details of our implementation of the 
delta field are relegated to the Appendix.) We discuss the power counting in Sec. IIIII In Sec. IIVI 
we show how several key ingredients of the ttN amplitude are calculated. The calculation and 
renormalization of the EFT amplitude are carried out in Sec. [V] We then fit the P-wave phase 
shifts in Sec. IVII A summary of our results and a conclusion are offered in Sec. IVIIl 

II. EFFECTIVE LAGRANGIAN 



In this section we briefly review how the effective Lagrangian is constructed. Much of this has 



already been discussed in the literature (see, e.g., Refs. [3, |6j, |7|, [L5J] ) , our review here serving mainly 
to establish the notation. Some details about the delta can be found in the Appendix. 
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A. Fields and symmetries 



The effective Lagrangian, with hadronic degrees of freedom, is expected to exhibit the symme- 
tries of QCD, which include Lorentz invariance, (approximate) parity and time-reversal invariance, 
color gauge symmetry, baryon-number conservation, and approximate SU(2)l x SU(2)r chiral 
symmetry. Although it is possible to extend the theory to incorporate violations of parity and 
time-reversal, we neglect them here since they are small relative to the accuracy pursued by us. 
Generalization to SU(3)l x SU(3)r is also possible. 

Color gauge invariance is trivially satisfied in the EFT because hadrons are color singlets. 
Baryon number conservation requires baryons to appear in bilinears. In the kinematic region 
where the EFT holds, external momenta are much smaller than the nucleon mass, Q <C m^, 
and thus Lorentz invariance can be fulfilled perturbatively in powers of Q/tun- There are two 
approaches in the literature to build the corresponding effective Lagrangian. One approach is 
to write a relativistic Lagrangian and then to derive the Q/tun expansion b y in tegrating out 



14 



20( | . one uses 



high-energy components using the path integral [30]. In the case of the delta 
a Rarita-Schwinger field. The so-called "off-shell" parameters that control the spurious spin-1/2 
sectors of the Rarita-Schwinger field are interpreted as choices of "gauge" . By building a gauge- 
invariant Lagrangian and choosing a certain gauge in Feynman rules, the spurious spin-1/2 sectors 
can be removed from the final result. 

It is not, however, inevitable to rely on the form of a Lagrangian outside the regime of validity 
of an EFT; only the symmetries should be important. Since in the region where ChPT is valid the 
nucleon and delta are always nonrelativistic, the other approach 
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151 ] starts from the nonrela- 



tivistic limit. This is accomplished by using heavy- fermion fields N for the nucleon and A for the 
delta, which are, respectively, two- and four-component spinors in spin (S) and isospin (I) spaces. 
Compared to the relativistic fields, the heavy- fermion fields have the common, inert, large mass 
removed, and contain only the destruction of particles. Not only does the heavy-baryon formalism 
keep clear track of the small expansion parameter, Q/mN, but it also is convenient as a framework 
to write the most general effective Lagrangian, where the only baryon degrees of freedom are those 
representing forward propagation. 

As usual, we introduce the spin operators S^ s \ normalized so that 

and in isospace t"', normalized the same way. We write = a/2 and t^) = T /2 in terms 
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of Pauli matrices. One also needs transition matrices that have proper Clebsch-Gordan (CG) 
coefficients embedded. We define 2x4 matrices Si in spin space such that their matrix elements 
between a nucleon state with a spin z-component a and a delta state with a spin z-component s 
are 



(i^H^;^) 



(3) 



where we use the notation of Ref. 



31]. It is not difficult to show that the bilinear N^SA is a 



three-vector, thanks to the Wigner-Eckart theorem. We impose the normalization condition 

1 



Si Sj 



(2Sij - ieij k a k ) 



(4) 



In addition, there is a spin-2 bilinear, N'QijA, which is a symmetric, traceless three-tensor, with 



1 13 1 

( 2 2 ;°"^'l 2 2 ; 2 ^ = 2 ^° iSi + CTjSi ^ s 



(5) 



Similar transition matrices, T and H a b, can be defined in isospace such that N^TA and N^E^A 
are an isovector and an isotensor, respectively. 

In order to build the effective Lagrangian, one first enumerates all the rotation-invariant opera- 
tors and then uses a Lorentz transformation rule perturbative in slow-velocity boosts to constrain 
the coefficients of those operators. However, it is not yet clear how one can go beyond order Q/win 
with this "bottom-up" approach. A separate paper by one of us [321] addresses this issue: the sys- 
tematic construction of a deltaful, heavy-baryon chiral Lagran gian with arbitrarily high relativistic 



321 ] relevant for the current paper 



corrections. The Appendix summarizes the aspects of Ref. 
In essence, we implement baryons via the Foldy-Wouthuysen representation [33] of the Poincare 
group, whose boost generators can be readily expanded in powers of Q/rriN- 

Chiral symmetry is more complicated to implement because it is spontaneously broken by the 
QCD vacuum to its isospin subgroup, and it is thus nonlinearly realized in terms of pion and 



baryon fie 
field b, 151 



ds jj, 



34 



351 ]. Here we use stereographic coordinates 7r to represent the isovector pion 



341 ] . for which we define a covariant derivative 



2U 



with f n ~ 92 MeV the pion decay constant and 



D = l + 



71" 



2 ' 



(6) 



(7) 



The nonlinear realization maps 



axial chiral-rotations of N, A, and into 7r-dependent (hence 



'local") isospin rotations [s,[34|,[35|]. Since such local isospin rotations do not commute with normal 



derivatives, one also needs covariant derivatives for N, A, and D,,. 



or a generic chiral-covariant 



field with isospin /, ip^\ its covariant derivative is defined as 

%^^(d, + tW.E^^ , (8) 

where 

E„ = i^-xD^. (9) 
For an isovector with Cartesian indices like D u , it is conventional to write the covariant derivative 



as 



%D U = 3 li D v + iE^ X D u . (10) 

Any isospin- invariant operator made of N, A, and their covariant derivatives will automatically 
be SU(2) L x SU(2) R invariant. 

Explicit chiral-symmetry breaking induced by the quark masses can easily be incorporated in 



the effective Lagrangian. T 



have a structure as follows 



rose oper 
3, Q,Q 



rose operators, denoted by that are proportional to m u ± will 



<f + = -D-'^-^ + D-^l-^j^, (11) 

.-=(*-^.,-) + ^-, ,12) 

where the quantities rj^ and r]f are built of covariant fields and are isovector and isoscalar, respec- 
tively. To preserve parity, it is easy to show that rj + (r]~) and rf^ (77^) are pseudoscalar (scalar). 
It can be shown [15[ that isospin is an accidental symmetry, in the sense that it only appears in 
the subleading effective Lagrangian. As a first study, we focus here on the isospin-invariant part 
of the ttN amplitude. 

Electromagnetic interactions can be easily incorporated in the Lagrangian by adding the re- 
quirement of U{\) gauge invariance. This is accomplished by turning all derivatives in existing 
interactions into gauge-covariant derivatives, and by adding additional gauge- invariant interactions 
with the electromagnetic field strength. Here for simplicity we neglect electromagnetic interactions. 



B. Effective Lagrangian 



Since the effective Lagrangian has an infinite number of interactions, one needs a scheme to 
organize all its operators. It is convenient to order the Lagrangian terms according to the so-called 

7 



chiral index v 



Q. 



f 



u = d + m + n 5 + ^- 2>0, 



(13) 



where d, m, ns and / are the numbers of derivatives, powers of m n , powers of 5 and fermion fields, 
respectively. The lowest value of the index is a consequence of the pattern of chiral-symmetry 
breaking in QCD. 

In constructing the Lagrangian, we use integration by parts and field redefinitions to remove 
time derivatives on baryon fields except for the kinetic terms. The Lagrangian terms with the two 
lowest indices are given by lfj 

2 

m 



2 ^tt 2 + NUSqN + g A N*raN • D 



£(0) = 2f 2 D 2 

+A* (i% - 5) A + 4^A¥f )£(i)A • .D + h A (VtSA + H.c^j ■ D + (14) 



and 



£ (i) = N 1 



m 2 



iV + 2c^ir 2 N^N + At 
2mjv D 



2rriN 



cm: 



A + 2^ 



m. 



^7T 2 A"tA 



'.N^TS ■ @A + ff.c. ) • D + 



(15) 



while the next-higher index yields 
Ami o 



£(2) 



7T 



5 Atil 2 A+ ^ 



2m 2 N 



2m 2 N 



Si m N^TS • ^A + W^TSi@ m A - 2ie ij iN^T£l im @ j A) + H.c. 



2D 2 

+ J}A_ 

8m 2 N 

+d ^ (l - ^ U iV f TSA + F.c.) • £ + 



( N^TS9 2 A - N^TS ■ ®§a\ + H.c. 



D 



(16) 



Here, g A (g A ) is the u = axial- vector coupling of the nucleon (delta) and h A (d) is the v = 
(^ = 2) 7riVA coupling. These low-energy constants (LECs) are expected to be of 0(1/Mq CD ) 
but are not determined by chiral symmetry. We define the phases of pion and delta fields so that 
g A > and h A > 0. The Am 2 term provides a correction to the pion mass that is proportional 
to the square of the average quark mass (it is related to the 1% term in Ref. Q]); in the following, 
in order to simplify formulas, we absorb its contribution in m 2 . Likewise, the nucleon and delta 
masses receive at v = 1 contributions that are linear in the average quark mass, the respective 
sigma terms (see, e.g., Ref. 36]) denoted here by c and c A . With our choice of heavy- nucleon field 



we have already absorbed cm 2 in the nucleon mass rajy. Again, for simplicity, in the following we 
absorb the remaining mass contribution, (c A — c)m 2 , in 5. The remaining pion-delta interaction, 



together with a number of other interactions not shown, contributes to the order we work below 
only to a further renormalization of 5. The interactions associated with the pion and nucleon 
mass corrections only contribute to our reaction at higher order. Note that "• • • " refer to terms 
that do not appear explicitly in ttN scattering to the order concerned in this paper. Higher-index 
Lagrangians can be constructed with more derivatives, but will also only contribute at higher 
orders. 

Different versions of the heayy-baryon effective Lagrangian that are deduced from a relativistic 
formalism are given in Refs. [l^.[l8|. Their Lagrangians, in our notation, both have an independent 
v = 1 irNA coupling (denoted by 63 and 63 + respectively in Refs. and [3]). This is 
because redundancy due to baryonic equations of motion is only removed at the relativistic level 
m Refs. QQ, and further minimization of the number of interactions due to the heavy-baryon 
equations of motion is not considered there. The d term in Eq. (1161) is equivalent to the combination 

n 

of couplings —2/4 + /s in Ref. [181 ]. 



III. POWER COUNTING 



We now turn to the ordering of contributions to physical processes. For definiteness we take 
ttN elastic scattering, although the power counting is the same for other one-nucleon reactions 
where the external CM energy can be dialed to near the delta-nucleon mass splitting. We consider 
throughout the case Q ~ ~ 5 <C Mqcd- 



A. Away from the resonance 

We first consider CM energies much below the delta peak. In this case power counting is 
standard [5- 7] with the simple generalization that 5 counts as Q. The contribution of a diagram 
with A nucleons (here A = 1), L loops, and Vi vertices with chiral index vi is proportional to Q p , 
with 



P 



2 - A + 2L + > 2 ~ A ■ (17) 



The contributions with minimum p form the leading order (LO), the next contributions are referred 
to as next-to- leading order (NLO), and so on. 

The power counting can also be applied to sub-diagrams if A is generalized to count any fermion 
line that is unattached on one side. Examples, which will be important later, are the following: 
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FIG. 1: The LO pion self-energy, . A dashed line represents a free pion propagator. The unmarked 
vertex has v — and the twice-circled vertex has v = 2. 




+ 



FIG. 2: The LO nucleon self-energy, H N . A single (double) line represents a free nucleon (delta) propagator. 
Vanishing diagrams are not shown. 

(i) The LO contribution to the pion self-energy E^ = 0(Q 4 /Mq CD ), see FIG.[H two powers down 
compared to the inverse of the free pion propagator. From the LO pion self-energy we obtain a 
correction to the pion-field renormalization constant, Z^ = 0(Q 2 /Mq CD ). 

(ii) The LO contribution to the nucleon self-energy E^ = 0(Q 3 /Mq CD ), see FIG. [21 two powers 
down from the nucleon kinetic energy. From the LO nucleon self-energy we obtain a correction to 
the nucleon- field renormalization constant, Zjp = 0(Q 2 /.Mq CD ). 

(Hi) The LO contribution to the delta self-energy E^ = 0(Q 3 /Mq CD ), see FIG. O two powers 
down from the delta kinetic energy and delta-nucleon mass difference. From the LO delta self- 
energy we can define a correction to the delta-field renormalization constant, = 0(Q 2 /Mq CD ). 

(iv) The NLO contribution to the delta self-energy E^ = C(Q 4 /M^ CD ), see FIG. H 

(v) The NNLO contribution to the delta self-energy E^ = C(Q 5 / m qcd)> see FIG.E1 

(vi) The NNLO contribution to the TriVA vertex form-factor vP = 0(Q 3 /M^ CD ), see FIG. El two 
powers down from the v = Ka vertex. 

This power counting holds for generic momenta, but it does not work equally well in every 
specific region of phase space. For example, close to threshold an incoming or outgoing pion has 
energy very close to m n and a three-momentum close to zero — it is, in other terms, nonrelativistic. 
In such cases, a spatial derivative on the pion field does not contribute the same as a time derivative, 
yet Eq. (|17l) does not discriminate between these derivatives. One can refine power counting for 




FIG. 3: The LO delta self-energy, Y, A . Vanishing diagrams are not shown. 
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FIG. 4: The NLO correction to the delta self-energy, S^- The once-circled vertices have v = 1. Vanishing 
diagrams are not shown. 
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FIG. 5: The NNLO correction to the delta self-energy, £V . The twice-circled vertices have i> = 2. Vanishing 
diagrams are not shown. 

this region if one wishes. 
B. Near the resonance 

It should thus be no surprise that the above power counting fails in the immediate neighborhood 
of the delta resonance. Consider the two contributions to ttN scattering in FIG. [7] at a CM 
energy E. Diagram (a) is proportional to 1/(12— S) and diagram (b) to T,^\e)/(E - 5) 2 , which 
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(2) 

FIG. 6: The NNLO correction to the irNA vertex, K ■ Vanishing diagrams are not shown. 




(a) (b) 
FIG. 7: Examples of one-A-reducible diagrams. 

are, respectively, 0(1/ Q) and 0(Q/Mq CD ) at a generic low energy. Although both E and 5 are 
separately of C(Q), they have opposite signs and can cancel. When that happens, diagrams (a) 
and (b) are enhanced above their standard power counting, and diagram (b) is enhanced more 
than (a) so that the self-energy might no longer be a small correction: a resummation is necessary. 



The situation here is completely analogous to other narrow resonances 26], where a "kinematic 
fine-tuning" requires a modification of power counting in the neighborhood of a resonance. 
In fact, within a window of size 

around the delta peak, the bare delta propagator is C(Mq CD /Q 3 ). Since is C(Q 3 /Mq CD ), one 
simultaneous insertion of and the bare delta propagator contributes 0(1), so the two diagrams 
in FIG. [7] become comparable. We should thus resum the geometric series of one-A-reducible delta 
propagators shown in FIG. Moreover, since within the window E = 5(1 + O (Q 2 / Mq CD )) , the 
energy dependence of T,^(E) is always two powers down from its value at E = 5: 

V%\E) = jl + O {jgrA | • (19) 

The resummation thus amounts to a dressed propagator 

S<£\E)= [E-S + E^iS)]' 1 , (20) 

which scales as Mq CD /Q 5 . This is an enhancement of two powers over the generic situation. To 
make a full amplitude, one needs to contract the dressed propagators with the ttNA vertices, for 
which we should also neglect the energy dependence: 

V^(E) = Vj n \5) L + O CfirA } • ( 21 ) 
12 
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FIG. 8: Dressed delta propagator at LO as a sum of insertions of £ 
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In contrast, in one-A-irreducible diagrams the delta propagators are far away from their pole 
and do not need to be dressed, continuing to scale as 1/Q. This is trivial in one-A-irreducible 
trees, e.g., the u-channel A-exchange diagram. Less trivial are the loop diagrams where there are 
integrations over the energy domain spanning the delta pole. In this case, pions in the loops carry 
at least m n of energy, and the delta will not go on-shell for ttN energies below 5 + m w . Even at 
this point, the delta will be recoiling and the kinetic energy needs to be taken into account before 
the self-energy. 

Take for example the one-loop diagrams with a delta that contribute to the LO nucleon (FIG. [2]) 
and delta (FIG. [3]) self-energies. Suppose that the external fermion and the internal pion have four- 
momentum p and I, respectively. Apart from the CG coefficients and coupling constants in front, 
the loop integral with the bare dressed propagator has the form 



One can always close the contour of integration over 1$ in the half-plane opposite to the half-plane 
where the pole of the delta propagator is. In this case, we pick the pole at Iq = oj — ie, where 



oj = Jl 2 + m 2 is the pion energy. The remaining integral over I has no singularities when po = 5. 
Alternatively, rewriting it as an integral over oj after carrying out the angular integrations, it starts 
at m n while the integrand has poles only at u = and oj = po — 5. The same argument holds for 
A = 1 diagrams with more loops. 

As a consequence, in one-delta-irreducible diagrams the standard ChPT power counting (|17p 
still applies; dressed delta propagators only need to be included in one-delta-reducible diagrams. 
We thus arrive at a new power counting for one-A-reducible diagrams within a narrow window 
around the delta peak, 



where «a is the number of dressed delta propagators. This is the non-electromagnetic version 
of p derived in a slightly different power counting in Ref. [^(J. We discuss the similarities and 
differences between the two approaches in Sec. IIIIDI 




(22) 





(23) 
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Notice that there is a larger region around the resonance, of size \E — S\ = 0(Q 2 /Mqcd), 
where the enhancement in the delta propagator is insufficient to compensate for an insertion of 
the self-energy. Although the power counting is slightly different than Eq. (|17fl . this case is still 
perturbative. For simplicity, we do not consider it separately from the generic situation away from 
the resonance, where \E — S\ = 0(Q). 

C. Sewing the two regions 

We can now weigh the diagrams contributing to 7riV scattering, putting together the two power 
countings (|17p in the off-the-pole region and (|23p in the pole region. These contributions generally 
scale as 0(Q P /M^ D f^), where 1//^ is due to the two pion external legs. Without causing con- 
fusion, we simply specify the size of the contribution of a 7riV scattering diagram by order . 
Found in FIG. [9] are diagrams up to order Q l . To this order, the kinematic fine-tuning to the 
delta simply promotes the diagrams on the right (labeled A, B, C) of the figure with respect to 
the standard assignment on the left (D). The diagrams at the bottom (E), which are one-delta 
irreducible, are not affected. 

It seems that the two different power-counting schemes, which are applicable in two different 
regions, would lead to an EFT amplitude in the form of a piecewise function in the energy. Even 
worse, separating these two regions is somewhat arbitrary. 

A piecewise EFT is actually unnecessary. Look at, for instance, the Q~ l order. Only the LO 
pole diagram (FIG. E|A)) contributes at this order. Since there is no diagram in the off-the-pole 
region, the EFT prediction should vanish when away from the pole. Enforcing the LO pole diagram 
in the off-the-pole region seems to make a "wrong" prediction. However, it is only wrong by an 
amount of order Q 1 , since a dressed delta propagator scales off-pole as a bare one (E — 5 dominates 
over ). This is of the same order as the tree in FIG. E[D). On the other hand, however, one 
already made an even larger error of order Q° by neglecting the NLO pole diagrams. To put it 
another way, extending the domain of the LO pole diagram is equivalent to shifting a subset of 
higher diagrams into the order Q . This does not disobey the original power counting as long as 
one does not claim a higher accuracy by doing so. In many perturbative quantum field theories, one 
runs the renormalization group to select an optimal renormalization scale in favor of a more rapid 
convergence of the perturbative series, which is also equivalent to a rearrangement of diagrams. 
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Away from the resonance Near the resonance 

\E-S\ = 0{Q) \E-8\ = <D(Q 3 /M>) 



"NT nnQ 

= in out? 


(A) 


Q° 

None 


(B) 




® © 


(D) 

Q 1 


® ® 


(C) 



FIG. 9: Contributions to 7riV scattering up to order Q 1 : (A) Q" 1 pole diagram; (B) Q° pole diagrams; 
(C) Q 1 pole diagrams; (D)fe(E) Q 1 tree diagrams, of which (E) apply to both regions. Power counting 
away from the resonance is the standard ChPT counting. The relation between the double line (bare delta 
propagator) and the thick solid line (dressed delta propagator) is given in Fig. [5J 

This sort of rearrangement, written symbolically as 

N a LO + a subset of ^ 1NPLO -> N a LO , (24) 

must be used with caution. First, the higher-order subset being added to N Q LO should be cutoff- 
independent by itself, in order to avoid introducing model dependence in the form of cutoff depen- 
dence. Second, one probably would not like to move undetermined LECs into lower orders because 
doing so weakens the predictive power of EFT. In our case, order- by-order renormalizability of the 
rearranged diagrams will be relatively simple to demonstrate because, in the resonance counting, 
energy variation in the delta self-energy and interaction vertices is treated perturbatively. We will 
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also see that the rearranged diagrams include in the off-pole region, within an error of order Q 2 , 
only delta parameters that appear in the diagrams of FIG. 0(D,E). 

Ensuring that our amplitude is correct also away from the resonance region does require some 
care. When the pion and nucleon in the initial and final states are on-shell, the momentum 
dependences in the ttNA external vertices translate into energy dependence. In the resonance 
counting, such energy dependence is subject to E = 6 + C(Q 3 /Mq CD ) so that the 7riVA external 
vertices in the LO and NLO pole diagrams are constant. This has two consequences, (i) Up to 
order Q° the P-wave amplitude in the delta channel does not vanish at threshold as it should. But, 
as we argued, applying the LO pole diagram near threshold is expected to make an error of order 
Q , which can be taken as "vanishing" in comparison with the dominant amplitude around the 
resonance, (ii) Since the energy expansion around 5 should not be enforced in FIG. E^D), simply 
continuing the LO pole diagram away from the resonance would not reproduce FIG. [9jD) due to 
the lack of energy dependence in the external vertices. To account for FIG. 0(D), the easiest way 
is to restore at order Q 1 the energy dependence of the external vertices of FIG. E5A). 

To summarize, with sufficient caution it is unnecessary to restrict the energy domain of the 
dressed pole diagrams in FIG. Of course, our power counting stresses the fact that amplitudes 
are much larger in the resonance region than at threshold. If one is interested solely in the region 
near threshold, where the delta enhancement is not relevant, one is better off with the standard 
ChPT power counting — or, for reactions where typical three-momenta are below the pion mass, 



with a "heavy-pion" EFT 



371 ] . Using our power counting in both threshold and resonance regions 



is only efficient if one aims for a unified description throughout the low-energy region. 



D. Other approaches 

The above power counting is not limited to ttN scattering. It applies to any reaction where 
one can dial the initial energy to bring the delta close to being on-shell. Of course, the need to 



account for the delta self-energy has been felt for a very long time, and simple tree-leve" 
with an added delta self-energy already provide at least qualitative descriptions of data 



models 



38|]. What 



EFT provides in addition is a way, consistent with QCD symmetries, to correct systematically for 
quantum effects and physics at short distances. We now compare our power counting to other 
approaches that considered non-perturbative effects of the delta in EFT. 

Reference resummed the LO delta self-energy in one-delta-reducible diagrams (as done 
here), and used the corresponding delta width as an independent empirical input. The scheme 
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is quite different from ours. First, it did not take into account the fact that the dressed delta 
propagator is enhanced by C(Mq CD /Q 2 ) in the resonance region, and hence adopted the standard 
ChPT power counting except for simply substituting the LO dressed delta propagator for the 
bare one. Therefore, the off-the-pole region was excessively emphasized with a two-order-higher 
accuracy than the resonance region. This discrepancy in power counting may be viewed as a 



rearrangement of the type (|24|) but Ref. JjJ enlisted many more LECs to achieve the same overall 
accuracy. Second, perhaps more importantly, Ref. [19( did not consider the corrections to the delta 

(2) 

self-energy. For instance, in our power counting the one-loop corrections to the irNA vertex, 
defined in FIG. [H contribute to the NNLO pole diagrams not only throu gh v ertices attached to 



191 ] . This would make 



(2) 

external legs but also in . The latter, however, were neglected in Ref. 
one, if using the scheme in Ref. [19], unable to recover Eq. (pQ) around the resonance, which is 
required by order-by-order unitarity. The difficulty of preserving unitarity around the resonance 
forced Ref. [3] to rely on certain prescriptions, referred to as S- or X-matrix method therein, to 
extract the phase shifts from the scattering amplitude. In particular, one of the prescriptions, the 
5-matrix method, leads to a discontinuity in the phase shifts. 

More sophisticated methods of unitarization based on EFT exist, for example Refs. 



24 



25|. In 



this case, if a sufficiently high-order kernel is fully iterated in an analogous fashion to what is done 



in the two-nucleon system 7J, 
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231 ] . corrections to the leading delta self-energy are accounted for 



and good results can be obtained even beyond the delta resonance. However, in a system with two 
heavy particles of reduced mass \x such a unitarization is justified by an infrared enhancement of 
0(47r /x/Q) over the standard ChPT counting [7]. While this might apply to the strangeness sector 



of ChPT, the strict EFT rationale for resummation in pion-nuc 



region must be rooted in arguments such as those in Refs. 



eon scattering in the resonance 



2CH22J and in the present manuscript. 



Contrary to Ref. 19|, the resummations in Refs. 2CH22f| are based on a power counting very similar 



to ours. Both in Refs. 



2014221] and here, it is recognized that the immediate vicinity of the resonance 



requires a different power counting than the one in standard ChPT. Moreover, in both approaches 
there is an attempt to smoothly bridge the resonance region with lower energies. 

Our 
in Ref. 



power counting differs from the one proposed in Ref. [20], however, in that 5 was assumed 
201 ] to be much larger than m n , rn^/MQCD ~ (S/Mqcd) 2 - As a consequence, the relative 



importance of explicit chiral-symmetry-breaking terms is reduced. Strictly speaking, in the light of 
the hierarchy m n -C 5, one would have to neglect the m 2 in the pion propagator in £ *^ and other 
diagrams where the pion momentum is 0(8). This might lead to unpleasant infrared divergences 
in certain diagrams, which Ref. [3] avoids by not enforcing their power counting in the calculation 
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of diagrams. Of course this can be justified by a rearrangement in the fashion of Eq. (|24p . and if 
this is done it blurs the difference between the two power countings somewhat. 

Clearly, as emphasized by the authors, the power counting of Ref. [3] is well-suited to study 
the regime of smaller quark masses, where m n decreases but 5 is basically unchanged. In contrast, 
it works less well as the number of colors N c increases, when 5 decreases. In any case, in the real 
world the scales are not clearly separated, 5 being larger than m n only by a factor of ~ 2, and the 
interesting limits — > and N c — > oo can always be studied separately afterwards. 

Thus, if one is going to generically treat the pion mass as comparable to momenta in loops 
around the resonance, we find it simpler not to emphasize such a factor of 2, given that there 
are lots of other similar factors floating around. We simply count 5 as comparable to the m w , 
?ti, t /Mqcd ~ 3/Mqcd, as it has been done before away from the resonance [13Hl8|. 



This discussion about the best way to power-count explicit-symmetry-breaking terms should not 
obscure — as it apparently has — the important fact that it is the kinematic fine-tuning to a narrow 
resonance (that is, one that has a width smaller than its energy) that demands a resummation 



[20! . This is in lad quite a general requirement of an LTT lor shallow resonances 
nothing do to with explicit chiral-symmetry breaking. 
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IV. ONE-DELTA-IRREDUCIBLE INGREDIENTS 

As an example of our approach, we want to calculate the 7rN T matrix to 0{Q 1 ). The tree 
diagrams are straightforward, but the pole diagrams in FIG. [9] are more complicated. The latter 
have three common components: the field renormalization constants and Zjy, the nNA vertex 
function V n , and the delta self-energy Sa- All of these ingredients are made up of one-A-irreducible 
graphs, which can be expanded in powers of Q/Mqqd according to the standard ChPT power 
counting. In this section we investigate each of the ingredients in turn, up to relative 0(Q 2 /Mq CD ). 

Note that in general these quantities are cutoff-dependent. If dimensional regularization (DR) 
is used, cutoff dependences in the loops appear as 1/(D — 4) poles in the complex dimensionality 
(D) plane near D = 4. Of course, the scattering amplitude should be independent of any cutoff. 
Therefore, the cutoff dependences must be absorbed by suitable counterterms. 
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A. Field renormalization constants 



Field renormalization constants Z^ and Zn are the residues of the corresponding two-point 
Green's function (or fully dressed propagator). 

The renormalized pion mass is the magnitude of the momentum of the pole of the dressed pion 
propagator, which we denote by m n . At the pole, the pion four-momentum p obeys po = &(\p\), 
the energy of an on-shell pion: 



uj(\p\) = yjp 2 + m 2 . 



(25) 



The pion self-energy, T,- K (p 2 ), is such that E 7r (m^.) = 0. The pion field renormalization constant 
Z v is related to T, w (p 2 ) by 



1 + 



_d_ 

dp 2 



(26) 



Expanding in powers of Q/Mqcd according to the standard ChPT power counting, to NNLO, 

4°) = 1 , (27) 



zp = 0, 



(2) _A S (0)( 2 ) 

dp 2 n U ' 



(28) 
(29) 



where t!§' is the LO pion self-energy shown in FIG. [TJ 

Similarly, we renormalize the nucleon mass in such a way that the pole of the dressed propagator 
of a nucleon of four-momentum p is given by po = E^(\p\), where 

(30) 



e n(\p\) = \ P 2 + m 



N 



m N 



is the heavy-baryon energy of an on-shell nucleon. The nucleon self-energy Y,n(p) satisfies 
T,n(En(\p\),p) = 0, and Z^ is related to Sat(p) by 



Z 



N 



1 + ^-E N ( P ) 

opo 



Po=E N (\p\) 



To NNLO, 



z (o) 
Z 



N 
(1) 



z 



(2) 
JV 



1 , 

0, 



dp 



(31) 

(32) 
(33) 
(34) 



Po=0 



where is the LO nucleon self-energy shown in FIG. El 

(2) (2) 

The field renormalization constants are not directly observable. Not surprisingly, Z N and Z\ 
are cutoff-dependent. 
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B. A propagator 



The loops appearing in the delta self-energy, E A (p), have both real and imaginary components. 
The inverse delta propagator is written as 

= Po - E A (\p\) + ReS A (p) + l 7 (p) , (35) 

where 



Ea(\P\) = ^P 2 + m\-m N (36) 
is the heavy-baryon energy of the delta, and 

7 (p) = 2ImE A (p) . (37) 

We choose to renormalize m A so that 

ReS A (£; A (|pl) ! p) = 0. (38) 

The regulator dependence that arises in the real part of the self-energy E A (p) at the on-shell point, 
Ej\{\p\), is thus absorbed in 5. 

Although the delta need not be an asymptotic state at least in low-energy irN elastic scattering, 
it is useful to introduce the field renormalization constant Z A , which we define as 

d_ 

Since the energy dependence of E A (p) only appears at NNLO, the leading few orders of Z A are 
given by 

4 0) = 1 . (40) 
Z« = , (41) 



Z A 1 = 1 + — [ReS A (p)] 



(39) 

po=e a (\p\) 



Z 



A 

(2) _ d 



ReE^(p) 



(42) 

Po=S 



A 0Po 

where Ea^ is the LO delta self-energy given in FIG.O We stress that Z A is merely an intermediate 
quantity in our calculation and the definition above is not unique by any means. As we will see in 
Sec. IV Bl the amplitude in the end does not depend explicitly on Z A . 
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C. 7riVA vertex 



We define the 7T./VA vertex function, V n , as the sum of all amputated A — > irN subdiagrams that 
have an incoming delta carrying four-momentum (po , p) and an outgoing pion with four- momentum 
(fcoj k) and isospin a. The incoming and outgoing particles are not necessarily on-shell. Rotational 
and isospin invariance require that V^a be a 2 x 4 matrices in spin and isospin space of the form 



S ■ kF + S ■ pG + €iji£li m kjpi (k m H + p m Q) , (43) 



where F, G, H, and Q are three-scalar form factors that depend on po, ko, p 2 , p ■ k, and k 2 . For 
simplicity, we study only the special case where the incoming delta sits in its own CM frame, i.e., 
p = 0. The vertex function is simplified as 

V™ = T a S-kF(p ,ko,k 2 ) , (44) 

with k = \k\. 

The function F at LO and NLO, F^ and F^ respectively, is constant and can be read off 
directly from C^> QID and £W (USD- At NNLO, F^(p ,k ,k 2 ) consists of v = 2 interactions 
and one-loop diagrams, shown in FIG. [6) F^(po,ko,k 2 ) has cutoff dependences proportional to 
Pq, k 2 ,, poko, 5po, 5ko, 5 2 , and m 2 . Naively, one could introduce v = 2 counterterms that are 
associated with time derivatives on A and/or 7r, in order to absorb divergences associated with 
Pq and ko. However, those counterterms can always be removed using equations of motion or 
field redefinitions, that is, they are redundant parameters, which S-matrix elements do not depend 
upon. Fortunately, the 7T./VA vertex function need not be cutoff-independent. Since there are other 
cutoff dependences floating around, such as those in the field renormalization constants and self- 
energies, what matters is that the combined cutoff dependence cancels out when the ingredients 
are put together in the EFT amplitude. 

As we will see in Sec. |V] a cutoff-independent form factor, Fr, will appear in the amplitude, 



Fr = VZWZNVZAF(5,oj(k s ),ki) , (45) 
where k$ is the momentum of the pion when all particles are on-shell, 



i 



It is worth noting that Fr is independent of energy by definition. Since the definition of Za is not 
unique interpreting T a S ■ k Fr as the CM amplitude for delta decay into irN is questionable. 

21 



We stress that Fr is just another intermediate quantity that will be used later to assemble the ttN 
scattering amplitude. 

The LO and NLO of Fr are the same as those of F: 



R ~ ~ 2f n 
F^ = = . (48) 



F (o, = F(0 ) = {A1) 

"JIT 



At NNLO, the pion loops bring dependences on 5 and kg. To the first approximation, 5 = Lu(k$) + 
0(kg/m]\r). In addition, F R 2 ^ receives nontrivial one-loop corrections, which are evaluated on-shell 
by definition, and contributions of the field renormalization constants, 

Ff = \ (z« + zf + Z%) F® + F^ (.(k s )Mk S ), kl) = (x + iX) , (49) 



where 



and 



with 



x= 5 



(4vr/ OT 



hA,R kg 



dR i — j§ + ReG(m n /ks) 



(50) 



fc2 

5 Im0(m„/fc) , (51) 



\2" 



= ^(l + ^)" l {-^(^-^f)^ 3 + 2 7 ri( 3 i + ^) 



+ 



5a-^^a, R (13 + 15x 2 ) + ^ 2 



Vl + a; 2 + 1 

There are two types of divergences arising from the one-loop corrections and field renormalization 
constants: one is proportional to 5 2 and the other m 2 . We can use the bare Ha to absorb the S 2 
divergence. In constructing (fTB]) . we already used the delta equation of motion to turn the 
operator -(N^TS&qA + H.c.) • 3 into 5 2 (N^TSA + H.c), which is subsequently absorbed into 
the Ha operator. Therefore, it is appropriate to combine the 5 2 divergence with the bare Ha, since 
there already is a 5 2 piece in Ha- A similar argument holds for the d operator, which is proportional 
to m 2 and can be renormalized by the m 2 divergence. Equations ([4T|) . (fl"9"]). and (f5Q|) should be 
viewed as the definitions of the renormalized coupling constants h,A,R and d,R. To make expressions 
compact, we will drop the subscripts "r" on Ha,r and cIr in the rest of this paper. 

Note that the logarithm in the function Q{x) blows up as x — > 0, implying that an infrared 
divergence would arise if one treated as a higher-order effect, as in the power counting of 

Ref. [20j]. Although this infrared divergence is not a fundamental difficulty, it is convenient to 
avoid it by considering m n ~ S. 
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D. A self-energy 



In this paper we employ the delta rest frame, where the on-shell delta energy is 

E A (0) = 5 . (53) 

As we discussed, the delta self-energy produces a large effect only inside the resonance window, 
\E — 5\ = C(Q 3 /Mq CD ). We will enforce this relation when expanding T, A (E). The restriction 
on the energy domain changes the order- index of T, A (E), since the energy dependence of T, A (E) 
around E = S is always two powers smaller than T, A (5), Eq. CE]). We denote by S^, S^, . . . the 
expansion of T, A (E) within the resonance window. 

The diagrams contributing to £a at LO, NLO, and NNLO can be found respectively in FIG.[3l 
FIG. HI and FIG. El Detailed calculations of the delta self-energy can be found in the literature, 



for example Refs. 



13 



39(]. Using Eqs. ([37}, ([38]), gSJ, and (|20D, the three lowest orders of T* A (E) 



can then be written as 



$2(6) = 












sf\E) 



(54) 
(55) 



+ -^ 2 HE), (56) 



where 



f 2)(£) = 7 (2) (ff ) + 4 2 ) 7 (0)(5) + (E- 5)7 (0) V) • (57) 

(2) - — ("21 

Due to the presence of Z\ , is cutoff-dependent. We should appreciate the fact that the energy 
dependence of ~ReT>^(E) is not present until NNLO, thanks to enforcing \E — 5\ = 0(Q 3 /Mq CD ). 
Had it shown up already in S^? , an otherwise redundant operator iA^^A would have been 
necessary in order to absorb the divergence in T,^(E) that is proportional to E s . 

The imaginary part of the delta self-energy, j(E), can be most conveniently evaluated by cutting 
the intermediate states that could be on-shell with an injected CM energy E and then replacing 
those propagators with the Dirac delta functions that enforce energy-momentum conservation. 
This is of course equivalent to applying the optical theorem. 

In LO and NLO delta self-energy diagrams, the only potentially on-shell intermediate state is 
a pion and a nucleon. The contribution of such an intermediate state to 'y(E) can generally be 
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written as 



1-kN 



(E) = Z N Z n J2 J 



^ ' V™\E)V™(E)x2ir5(E-E N (l)-u(l)) 



(2tt) 3 2w(I) 



A' 8 



Z N Z n Af(k)— \F (E,uj(k),k 
Git 1 v 



2\ |2 



where, in the last line, k satisfies 



and the pre-factor is 



E = uj{k)+E N {k) 



_ E N (k) + m N 
E + rriN 



For an on-shell delta, use of Eq. ([4"5]) yields 

7ffJV («5)=^)^lM 
Substituting the LO expressions for F (Eq. (|i7]> ). ^ (Eq. ([271)) and Zjv (Eq. ([32D ). 

7 (0)(£) = jL_X(k)k 3 . 



(58) 



(59) 



(60) 



(61) 



24tt/ ; 



(62) 



Using in J\f(k$) kg the exact kinematic relation between and 5, Eq. (06]), 

(*» - ml)i [1 + */nw 4- - »4)/(*»»n 1 1+V "^ + '^;"' )/2 '"" 



7^(5)= ^ 



24*/, 



a relation known from isobar models 
Z n (Eq. d28j)), and Z N (Eq. 



38] 



(63) 

Analogously, from the NLO expressions for F (Eq. (j48[> ). 
(1) (<5) = 0. (64) 



7 



The strict heavy-baryon expansion of Af(kg) kg gives 



7 (O)0) +7 W ($)+... 



Ml 



^2 m 2N 5 



247T/I 



55 , 42<5 2 - 7m 2 7£ (225 2 - 7m 2 ) 



2l7lN 



8m 2 N 



(65) 



The first term is the well-known heavy-baryon limit [13]. With 5 ~ 300 MeV, the first, Galilean 
correction is ~ —80% due to the relatively large numerical factor in front of 5/mjv- This expansion, 
nonetheless, still converges as long as 5 < mj\r- More importantly, the slowness of this expansion 
is not reflected in the EFT expansion of the amplitude. As we will show later, to NLO the ttN 
amplitude only depends on S and 7^°- ) + 7*- 1 -*, i.e., it is a simple Breit-Wigner formula. Enforcing 
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the heavy-baryon expansion in N(k$) kg does not result in an amplitude with a different functional 
form, but it does lead to a different value of Ha- This only means that the value of Ha depends, 
because of the slow convergence, on how Q/rriM corrections are treated. In addition, there are 
higher-order corrections to 7 contributed by the NNLO 7riVA vertex function, which includes 
one- loop corrections and an undetermined LEC (d). Since there is no a priori evidence that these 
undetermined corrections are smaller than those that are proportional to l/m? N and higher, keeping 
Q/mjs! terms to all orders does not have any deep significance. 

In the following, we consider two cases. In one case, having explained the validity of the heavy- 
baryon formalism despite the large factor of —55/2mj^, we carry out the expansion of M(k)k s . We 
refer to this as the strict heavy-baryon expansion. In the other case, we do not expand J\f(k)k 3 , 
because (i) it is a convenient way to include the required terms up to l/m? N ; (ii) it allows a 
meaningful comparison of LECs with the literature, where a similar resummation is performed 
HQ; and (Hi) whenever desirable, a strict heavy-baryon expansion in Q/tjin can easily be 
worked out. We refer to this second case as the semi-resummation. 

When evaluating 7^ 2 ^((5), there are a few NNLO delta self-energy diagrams in which an interme- 
diate state of two pions and a nucleon (ttttN) could be on-shell if being cut on the middle-nucleon 
internal line. These diagrams are labeled (a) to (h) in FIG. To estimate the contribution of 
7T7riV to j^(5), we first notice that the phase space for such an intermediate state to go on-shell 
is so small that even pions are nonrelativistic, having three-momenta Q ~ y / m 7T (5 — 2m„) ~ 40 
MeV <C m n . We can use this fact to refine once again the standard ChPT power counting. To be 
definite, let us look at the two-loop diagram labeled (a) in FIG. In contrast with the generic 
situation, we should now replace Q for Q where pion three-momenta appear in the loops (vertices, 
propagators, integrals). Also, the energy of the nucleon is of 0(Q 2 /mjy) rather than Q. Overall, 
these changes in the power counting bring a suppression of roughly order (Q/m^)^ 7 ~ — 10 . 
This somewhat crude estimate is justified by a phase-shift analysis [28]. The unitarity of the S 
matrix suggests that the opening of the A -4 ttttN channel brings "inelasticities" in ttN scattering. 
However, Ref. 28] gives inelasticities only of order of magnitude 10~ 3 in the delta region. There- 
fore, though they are formally NNLO, the ttttN contributions are suppressed by the "accidental" 
closeness of 5 to 2m n . Numerically we can thus safely neglect ttttN contributions, 

IttttN^) oi , (66) 

and consider only ttN contributions. 
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With this approximation and Eq. (09]), ^ 2 \S) is given by 

^\5) = {2x-Z^)^\5), (67) 
where x is given in Eq. (foT)|) . Inserting Eq. (|6T|) into Eq. ([FT]) one finds 

7( 2 )(,B) = 2«7^(<5) + 7 (0) '(<5)(£ " *) , (68) 
completing the calculation of the delta self-energy to NNLO. 

V. irN SCATTERING AMPLITUDE 

In this section we put together the various one-A-irreducible ingredients calculated in the pre- 
vious section. We first review the kinematics in ttN scattering, before constructing the T matrix 
in the various channels. 

A. Kinematics 

In the CM frame, we denote the initial (final) pion momentum by k (k'), the initial (final) 
pion isospin index by a (a'), and the initial (final) nucleon spin z-component and isospin third- 
component by a (a') and r (r'), respectively. The CM energy, denoted by E, is given in terms of 
the pion (nucleon) energy oj(k) (Ejy(k)) by Eq. (JSSDi with k = \k\. 

The T matrix is related to the S matrix by 

S = 1 + %T . (69) 
Asymptotic pion states are normalized so that 

(iv,k' a'\n,ka) = 2ui(\k\) (2tt) 3 5^(k - k') 5 a i a , (70) 

while for nucleon states, 

(N, k' a' t'\N, kcrr) = (2vr) 3 5 {3) (k - k') 5 a ! a 5 T ! T . (71) 

The scattering amplitudes, A. a > a , are the elements of the T matrix between the asymptotic pion 
and nucleon states and can be written as 2 x 2 matrices in nucleon spin and isospin indices, 

(A a , a ) al(riT , T (k',k) = (k'a'a'r'\T\kaaT) . (72) 
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We normalize the spin-orbital projector Fji for a total angular momentum j and an orbital 
angular momentum I so that 



dn y ,F jH ,(k', k'^jiik", k) = 5 jfj 8 v iPji(kf, k) , (73) 
with dVty, the area element on a unit three-sphere spanned by k" . The Pj/s with lowest js and Is 



arc 



°i {k',k) - 4yr , 



1 l>/ 



4tt 



k'-k + 



(k' x k*J ■ <j 



Ps, (&', fc) = — fojfe' ■ k — i (k 1 x k) ■ a 
a 1 47T L V / 

The isospin projector It for a total isospin t is normalized so that 

} j It'(a , ,c)I t (c,a) = 5 t 'th(a,a) . 

c 

There are only two different I^s in ttN scattering: 

Il(a',a) = \{5 a 'a + ita'acTc) , 
2 3 

Ilia', a) = -(25 a > a - i€ a > ac T c ) . 

2 3 



(74) 
(75) 
(76) 

(77) 

(78) 
(79) 



The normalization factor between the angular-momentum eigenstates and the asymptotic states 
defined in Eqs. ([70]) and ([7T]) can be found in many textbooks (see, e.g., Ref. Q]). Without showing 
the tedious details, we simply state that A a ' a is related to the phase shifts, 9ju{E), as follows: 

iAafatf, k) = -^^y F Ak', k)l t {a', a) {exp [2i9 jlt (E)} - 1} , (80) 



jit 



where J\f(k) is defined in Eq. (|60p . The partial- wave T-matrix elements are expressed in terms of 
Ojit(E) as 



T jlt {E) = -i{exp [2i9 jlt (E)) - 1} . 



(81) 



In the following we will use a more conventional notation for a specific partial wave: ^2t,2j- For 
example, -P13 refers to the I = 1 (P wave), i = 1/2, and j = 3/2. 



B. T matrix 

Now we collect all the pieces from Sec. HVI to build the irN scattering amplitude in the various 
waves, Eq. ([8T]) . Here the exact relation between E and k (Eq. ([59]) ) is assumed, meaning that 
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certain trivial, kinematic k/m^ terms are resummed — what we refer to as semi-resummation. In 
the next section, Sec IV CI we specialize to the strict heavy-baryon expansion. 

At LO there is only a pole diagram, FIG. [9^A), which contributes only to the -P33 wave. 

From Eqs. (M and (Rg 



T-.LO 
1 P 3 3 



(0) 



1 + 



Q 



Mi 



QCD 



(82) 



E-8 + vyW(8)/2 

where j(°'[6) is given by Eq. (|63p . This is of the form (JTJ) with a resonance at .Er = 5, an energy- 
independent width r = 7^) (5), and no background, Tg = 0. The two independent parameters can 
be taken to be 5 and h^_. This is same result as in any isobar model with the simplest contribution 
to the width resummed 38]. 



In next order (Q°), there appear Q/rriN corrections to the pole diagram, FIG. [9](B), which 
contribute via corrections to the delta self-energy, Eq. (|55p . The NLO amplitude has the same 
form as LO, 



rpNLO 

J p 33 



7 (o)(£) + 7 W(6) 



(l)f 



l + O 



<3 2 



^QCD 



(83) 



^-5 + i[ 7 (0)( ( 5)+ 7 (i)(<5)] /2 

where 7^(5) vanishes, as given by Eq. (|64p . when we do not expand kinematic relations in powers 
of 6/rriN- 

The NNLO (Q 1 ) corrections are more complicated. The cutoff dependences in \ Zj^\ 
(2) 

and Zn make several NNLO pole diagrams (labeled from (1) to (4) in FIG.[UJC)) divergent. The 
remaining diagrams in FIG. EJC) involve only Q/m^ corrections. To see that the divergences of 
these diagrams in fact cancel each other, consider the sum of these pole diagrams, 



nPole(2) 

-P33 



7 (0) (*)S£ 



(0) 



^ (2) (co(k s )Mk S ),kl) -S^^(E) 



-yW (S)(E - 6) 
7 (o)(tf) 



(84) 



Using Eqs. (09]), (f56j) . and ([68]) . and defining 



T B (5) = 2X, 



(85) 



we find 



r pole( 2 ) = _ 7 (0) ((5)s (0) 



2^ A 7 7(0) (5) 



+ iT B {5) 



(86) 



where no cutoff dependence is present. 
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Now summing up all the pole diagrams up to order Q , 



ppole 

■P33 



rjiLO I rp. 
1 P33 +1 P 33 

- 7 (°)(«5)5. 



pole(2) 



(0) 
A 



2 5 * 7 (JB) + ^)(5)" 



+iT B (^) 



Q 3 



M QCD, 



(87) 



Within the stated error, we can resum the corrections, 



-ipolc 



7 (o)(«y)+^(2)(S) 



[l + ir B (5)] 



Q 3 



M Q CD/ 



Pm ' £- J + i [ 7 (°)(<J) +7(2)^)] /2 
The result so far relied on an expansion around the resonance, to be joined with a description 



of the off-pole region. As we discussed in Sec. IIII CI FIG. 0(D) can be accounted for off-pole by 
allowing the "full" energy dependence in the external legs in FIG. 0(A) . This eventually amounts 
to replacing in Eq. ([88]) 7 (0) (£) with ^(E) and, in addition, 7 (2) (£) with 7 (2) (<5) = 2x^(5) in 
order to avoid over-counting the energy dependence. Moreover, to the order we are working, we 
can instead replace j( 2 \6) with [(1 + x) 2 — l] 7 (°)(i2). We then arrive at 

Q 3 N 



rppole/ off-pole 



E-S + iT(E)/2 



[l + iT B (S)] 



l + O 



M QCD 



(89) 



where, using Eq. (|62|) . 



T(E) 



[h A (i + x)\ 



k 3 M(k) . 



(90) 



247T/2 

Equation (|89p resembles Eq. ([TJ) except for the absence of an additional background term Tb(5), 
which would sit outside the pole term. The one-A-irreducible trees at order Q 1 will provide the 
remaining piece expected from unitarity, as we will now show. 

In fact, in order to complete order Q 1 , we need to include the remaining trees, FIG. 0(E), which 
are one-A-irreducible. Their contributions to the P33 channel are found to be 



with 



Tb(E) 



T tree, lAI^^ 



k 3 M(k) (g\ , 1 h\ 



+ 



(91) 



(92) 



6tt/ 2 V£ 36£ + <V 
This reduces in leading-order in Qjm^ to Tb(<5) defined in Eq. (|85p when E = 5, which is exactly 
the missing piece expected from unitarity. This is certainly not a miracle as EFT is supposed to 
reproduce the unitarity and analyticity of the S matrix order by order. However, it does confirm 
the consistency of our power counting from a particular perspective. 
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Now one can sum up pole (|89p and tree (|9ip contributions to obtain the P33 partial-wave 
amplitude. Since the difference between E and 5 is higher order in the NNLO pole term 



rpNNLO 

± Pi3 



T(E) 



[l + iT B (E)]+T B (E) + 0\T}£ 



Q 3 



M QCD, 



(93) 



E-6 + iT(E)/2 

Again we recover a Breit-Wigner form ([1]), but now with an energy-dependent width, Eq. (|90p . 
and an energy-dependent background, Eq. (|92p . The NNLO amplitude involves four independent 
parameters, S, Ha, gA, and x. The EFT thus provides specific energy dependences for T and Tb 
through the two extra parameters that appear at NNLO, gA and x. Our amplitude is similar to 
the one from isobar models [38j, except that Ha in T{E) gets corrected by a factor of 1 + x. 



Equation ([93]) reduces to Eq. ([I]) with energy- independent T{$) and Tb(6) when \E — <5| is 
small enough. For that to happen, the contribution from Tg(<5) must overpower that of the energy 
dependence of T(E), which according to Eq. (|86|) requires the size of the energy window to be 
smaller than the half width. In terms of the scales in the problem, we can define a "Breit-Wigner 
window" \E — S\ ~ 0(Q 4 /Mq CD ) where Eq. dU) holds approximately. Outside this small window, 
the more general form (|93|) should be used. 

Other channels are easy to calculate from the one-A-irreducible tree diagrams in FIG.[9^E). For 
the remaining P-wave channels, 



thNNLO 

J Pl3 



t-hNNLO 
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12tt/| 9E + 5 
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QCD 



(94) 



The one-A-irreducible contributions, Eqs. (|92|) and (jMj) . are, of course, related to amplitudes 



found in the literature. They are part of the simplest tree-level delta-isobar model [38J, where the 
relations among the NNLO amplitudes in P13, P31, and P\\ were also noticed. These contributions 
are LO in both deltaless and deltaful EFTs when one uses the standard power counting. Ignoring 
the delta contribution, namely the terms proportional to h\, one reproduces the deltaless EFT 



results in Refs. [9 
EFT given in Refs. 



ldLlnc 



12 



17 



uding the delta contribution, our results agree with those in the deltaful 



The seagull diagram in FIG.[9jE) (the Weinberg- Tomozawa term) contributes to the £11 and S31 
channels. Since the delta does not contribute to these waves, our results reduce to the well-known 



ChPT expressions at standard LO 



. In the following, we focus on the P-waves. 



C. Strict heavy-baryon expansion 

We have partially resummed Q/m^ terms in the amplitudes shown above. More precisely, 
this semi-resummation consists of keeping the exact kinematic relations ()46|) between 5 and kg in 
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Eq. (j61]) . and (|59|) between k and £" in Eqs. (|90|) . (|92jh and (|94j) . If, instead, one enforces a strict 
1/to/v expansion, there are some small changes in the results. 

In the P33 channel, at LO and NLO the amplitude is still given by Eqs. ([82]) and ([83]) . except 
that 7 (°)(£) and 7 W(5) are given by the first two terms in Eq. (|65p . At NNLO, again the amplitude 
remains in the same form as Eq. (|93|) but with different expressions for the width, 



t(e) = 

and for the background, 
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6VT/2 



E 36 £ + 5 



(96) 



Like for Tb(E), in other P- waves we simply replace k with y E 2 
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- rpNNLO 
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QCD 



(97) 



12VT/2 

The difference between these expressions and the corresponding ones in the previous section is 
of higher order. Comparing the effects of the two sets of formulas gives a further estimate of the 
size of higher-order corrections. 



D. Scattering volumes 



Although we do not aim at a precise description of the threshold region, we can extract from 
our calculation the 5"-wave scattering lengths and the P-wave scattering volumes, which are related 
to the amplitudes at zero energy. 

Up to order Q 1 in the off-pole region, the scattering lengths are, of course, identical to the 



venerable current-algebra results 



a Pa: 



The scattering volumes, with the semi-resummation, are 
\ 2 



1 ' — 1 1 - ml/5 2 U + 



a Pl , 



a P .,_ 



(l + m w /m N ) 1 



24iTf 2 m^ 



y/2hj 



m T /5 



3g A / l + m^/5 



(1 + m n /m N ) 



-1 



(98) 
(99) 



In the strict expansion, the factor of (1 + m w /m^)^ 1 should be dropped. The extra suppression 
m^/5 of the delta contributions at threshold is evident, and when m^/5 — > the delta decouples, 
as it should. In the opposite limit, m^/5 —¥ 00, the delta contributions grow and the scattering 
volumes vanish (to order Q 1 ) if KaIqa = 3/\/2, as is the case in the large-iV c limit j^. (Of 
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course our results do not direct apply to this limit, where the delta and other degrees of freedom 
are degenerate with the nucleon and do not decay into a nucleon and a pion.) The real world is 
in-between these two limits: since m^/S ~ 1/2 and \iaI9a ~ 3/\/2 (as we are going to see), delta's 
contributions are neither terribly small nor completely opposite to the nucleon's. 



VI. P-WAVE PHASE SHIFTS 



To test our EFT result, we compare it with the phase-shift analysis (PSA) by the George 

n 

Washington (GW) group [28], which bridges over the delta resonance. With parameters extracted 
from the fit, we compare the resulting values of the scattering volumes with those obtained by a 



PSA that focuses on lower energies 



29] 



We first extract the phase shifts out of the EFT amplitude in such a way that unitarity is 
preserved perturbatively. For that, we use Eq. (|8ip . expanding both the phase shifts and the T 
matrix in powers of Q/Mqcd- denoting by a superscript ^ the corresponding power, 



exp 



(n) 



(100) 



Specifically, in the P33 channel at LO, 
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and at NNLO, 



aNLO 

^33 



oNNLO 

^33 



arccot 



V 0) (<5) 



S-E 



7 (o)($) +7 (i)($) 



arccot 



t 5-E 



+ 



T B (E) 



(101) 



(102) 



(103) 



In the other P channels, 



oNNLO 



gNNLO 



iNNLO 
Pu 



1 rpNNLO 
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(104) 



Several LECs enter these EFT results. A number of them can be determined from other 
processes, such as pion decay. We adopt the following values: m n = 139 MeV, ttin = 939 MeV, 
qa = 1-29, and f n = 92.4 MeV. The value for is obtained using the Goldberger-Treiman relation, 
gA = 4v / vr/ 7r //m 7r + , from the pion- nucleon coupling constant / determined very precisely in the 
Nijmegen PSA of two-nucleon data [4l|. Note that this includes the chiral-symmetry-breaking 
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corrections to the pion-nucleon vertex (the Goldberger-Treiman discrepancy), which only appear 
two orders higher than the highest order we are working at. We have verified that our results 
are not very sensitive to such small corrections, and that using the chiral-limit value does not 
affect any of our conclusions. In principle we could determine \gA/fn\ from our fit, but that would 
distract from our main objective, the description of the delta resonance. The LECs germane to 
delta-resonance physics are S, Jia/ fn, an d *c. These LECs are to be determined from low-energy 
reactions involving the delta, and we do so here for ttN scattering. 

Our strategy of fitting is to determine the free parameters, 5, Ha, and k from the P33 phase 
shifts around the delta peak and then predict the phase shifts at lower energies in all P waves. As 
it is clear from Eqs. (jlOip and (|102p . the LO and NLO -P33 phase shifts have the same functional 
dependence on E: LO and NLO only differ by Galilean corrections buried in the width. In the 
following figures we will display Eqs. (jlOip and (|103p by lines. We fit our curves to the results of 



the GW PSA 



281 ] . which we indicate by dots. The two (four) points around 9p. i3 = tt/2 used to 



determine the two (three) parameters at LO (NNLO) are explicitly marked. 

The difference between Eqs. (jlOip and (I103p is itself an estimate of the systematic theoretical 
error at LO. We estimate the error at NNLO by power counting the higher-order corrections 
that are neglected. As shown in Eq. (|93p . the systematic error in the P33 channel is of size 
0(T^Q 3 /M^ CD ), while in the other channels it is of O (T^° A Q / M QCD ) , see Eq. ([MD- To 
make the estimate more concrete, we take for Mqcd the scale associated with the lowest-lying 
baryon integrated out in the EFT, the Roper resonance of mass rriR = 1440 MeV. For simplicity 
we take the factor to be tjir — mjy in all channels; while this could be a conservative estimate of 
the error since Roper contributions are likely to be small in channels other than P\\, it might still 
be representative in the P33 channel, where the next resonance, the A(1600), lies above the delta 
a similar distance. Since P waves other than P33 are predictions that apply throughout the range 
of energies we consider, Q should be given by the pion momentum, so the errors in these channels 
can be generically estimated as 

AT Non _ A = ±T™ L _° k - . (105) 

m R - m N 

In the -P33 channel, on the other hand, we take Q as a measure of the deviation from the resonance, 
because we would like to make the estimated uncertainty vanish at E = 6 where we fit to the PSA 
inputs. We first estimate the errors of T(E) and Tb(E), then calculate the systematic error of the 
amplitude using Eq. (|93p . In doing so the unitarity condition is explicitly preserved. The errors 
of T(E) and Tg(E) should be proportional to E — 5. Near the resonance, where T(E) matters, 
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E - S ~ Q 3 / m qcd- The third-order correction to T(E) is of O(7 (0) (5) Q 3 /M^ CD ), which then 
could be estimated by 

AT(E) = ±T{E) E ~ 6 . (106) 
m R - m N 

The higher-order corrections to Tb(E) come from one-A-irreducible tree diagrams with one v = 1 
vertex, and hence are of 0{Tb(E) Q/Mqcb)- Since, like in other channels, these one-A-irreducible 
diagrams apply throughout the low-energy region, E — 5 ~ Q and 

AT B (E) = ±T B (E) E ~ 6 . (107) 
m R - m N 

The systematic-error bands of the NNLO results are shown in the following figures as shaded 
regions. 

The phase shifts for the strict heavy-baryon expansion are shown as dashed line for LO 
(Eq. (jlDip ) and solid line for NNLO (Eqs. (fTU5j) and ([TOIj) ) in FIGs. [TO] and EJ 1 The P 33 EFT 
results agree quite well with the GW PSA past the delta resonance, and then, as expected from a 
momentum expansion, start to deviate from the data. We see a reasonable convergence pattern, 
and the small error band helps us understand why the EFT works so well in that channel: there 
is not much room for the EFT to change at higher orders. The P13 and P31 data are not far from 
the equality predicted at NNLO, and the Ua coming out of the P33 fit is such that the EFT curve 
provides a good description of the average phase shifts, being slightly below the data in the P13 
channel and above in P31. In both cases the empirical phase shifts are within the error bands. In 



the three channels our fits are comparable to other unitarization methods [24 

nn 



25j. Well below the 

resonance energy, our results are comparable to others at 0(Q 1 ) [l7|, |l8j], but not as good as 0(Q 3 ) 



I18H . which contains further LECs. On the other hand, the discrepancy with data in P11 is out of 

... nn 

the error band and thus significant. A similar discrepancy was already seen at LO in Refs. [10|,ll8|, 
where the focus was energies below the resonance. The discrepancy could be considered a "small" 
effect compared to the dominant P33 amplitude. Indeed, Refs. fiol . found that higher-order 
corrections improve the EFT description in the P\\ channel, allowing a change in the sign of the 
derivative of the phase shifts at the price of fitting more LECs. Unitarization including the Roper 



24, 



25j also works well throughout the energy region we consider here. 
The LECs extracted from the P33 fit are given in TABLE HI where they are denoted "strict" 
to emphasize that they were obtained from an amplitude where higher orders in Qjm^ were not 



1 The LO curve in the P33 channel does not vanish at threshold, but, as explained in Sec. IIII CI the deviation from 
zero is higher-order compared to the resonance amplitude. Indeed, the P33 phase shift vanishes exactly at threshold 
starting at NNLO. 
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FIG. 10: P33 phase shifts (in degrees) as a function of Wcm (hi MeV), the CM energy including the nucleon 
mass. The EFT strict heavy-baryon expansion at LO (NNLO) is represented by the red dashed (black solid) 
line. The NLO curve coincides with LO. The light-blue band outlines the estimated systematic error of the 
NNLO curve. The green dots are the results of the GW phase-shift analysis Points marked by a red 
star (black square) are inputs for LO (NNLO). 



TABLE I: Low-energy constants extracted at LO, NLO, and NNLO from the fits using the strict heavy- 
baryon expansion (strict) and the partial resummation (semi). 





5 (MeV) 
LO NLO NNLO 


h A 

LO NLO NNLO 


X 

NNLO 


strict 


293 293 320 


1.98 4.21 2.85 


0.050 


semi 


293 293 305 


2.71 2.71 2.92 


0.058 



treated differently than higher orders in Q/Mqqt>- One can estimate the errors in the NNLO values 
as the variation in each LEC within which the NNLO P33 curve in FIG. [TPl roughly stays within 
the error band. This is of course not a rigorous statistical method; it only serves to indicate how 
confident we are about the fitted LEC values. This way we find <5/MeV, Ha, and k to be within 
~ ±4, ±0.30, and ±0.030, respectively, of the NNLO values in TABLE H 

The delta-nucleon mass splitting is related to the position of the delta pole, which can be found 
by seeking the root of 



S^{E) = E-5 + iT{E)/2 = . 



(108) 



It yields the values given in TABLE |TI] under the label "strict" , which agree fairly well with the 
values from the GW PSA 28| and from the Review of Particle Physics [1]. In addition to the 
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FIG. 11: Predicted phase shifts (in degrees) in the P13, P31, and P\\ channels as functions of Wcu (in 
MeV), the CM energy including the nucleon mass. LO and NLO vanish in these channels; NNLO EFT 
results in the strict heavy-baryon expansion are given by the black solid lines. The light-blue bands outline 
the estimated systematic errors of the NNLO curves. The green dots are the results of the GW phase-shift 
analysis [2^ , 



value of 5 in TABLE ffl we find at LO T(5) = 104 MeV and at NNLO T(S) = 246 MeV. The 

n 

Breit-Wigner values [1| <5bw ~ 1232MeV — tun and Tbw ~ 118 MeV are extracted from a fit of the 
form (pQ) to data around the resonance peak. Since it is not clear that this window coincides with 
the "Breit-Wigner window" where our results reduce to Eq. ([I]) , we cannot adopt 5bw and Tbw as 
values for 5 and T(5), even though the pole positions agree. When other reactions are considered 
within our approach, one should use the value of 5 determined above. The large value of T(5)/2 
at NNLO indicates that it is not a good approximation for the imaginary coordinate of the pole, 
because the derivatives of T(E) in Eq. (|108p are not small. Since our NNLO amplitude gives a 
good fit of the phase shifts, there is nothing intrinsically wrong with our T(E); what we see here 
is an example of the general argument given in Ref. Q] that F(6) is scheme-dependent. 

The NNLO value of Ha we found is consistent with other sources, for example: h A = 1.96 — 2.64 
Q, h A = 2.92 Q, and h A = 2.81 Q]. Moreover, we obtain UaIqa = 2.21, close to the large--/V c 
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TABLE II: Pole position of the delta resonance (in MeV) extracted at LO and NNLO from the fits using the 
strict heavy-baryon expansion (EFT-strict) and the partial resummation (EFT-semi) compared with values 
from the GW phase-shift analysis (PSA) [28( and the Review of Particle Physics (RPP) 1|. EFT results at 
NLO are identical to LO. 





LO 


NNLO 


EFT-strict 


1232 - 52i 


1211 - 50i 


EFT-semi 


1232 - 52i 


1211 - 48* 


PSA 


1211- 49. hi 


RPP 


« 1210-50i 



TABLE III: P-wave scattering volumes in units of m~ 3 (with the charged pion mass): EFT in a strict 
heavy-baryon expansion (EFT-strict) and with partial resummation of relativistic corrections (EFT-semi) 
compared with results from a partial- wave analysis (PSA) [291 ] . 





a P 3 3 


a Pl3 


a Psl 


«Pii 


EFT-strict 


0.20 


-0.034 


-0.034 


-0.13 


EFT-semi 


0.18 


-0.028 


-0.028 


-0.11 


PSA 


0.2100(20) 


-0.03159(67) 


-0.04176(80) 


-0.0799(16) 



0. 



ratio h^/gA = 3/\/2 [40J. The third parameter in the fit, x, only appears at NNLO, and has 
several LECs embedded in itself, see Eq. (|50p : d and cn£, which have not yet been pinned down 
at the order we consider here. The central value we obtain is very close to the naive estimate 
(5 2 — m^)/(47r/„-) 2 ~ 0.05, but from the estimated error we see that x is not determined precisely 
in our fit. 

Substituting the fitted values for h A and 5 at NNLO in Eqs. ([35]) and ([99j), we find the EFT 
predictions for P-wave scattering volumes labeled "strict" in TABLE IIIII The ratio hj^jgA over- 
compensates for the m n /5 suppression, so that the delta contributions are not negligible. The EFT 



predictions are compared to the values extracted from the low-energy PSA of Ref. |29j, which are 
consistent with earlier extractions. We see good agreement for P33 and P13, and P31 is not too far 
off. The 1/4 ratio between Pn and P31, predicted in Eq. (|99p . is not respected well in the real 
world. This is a reflection of the discrepancy between the EFT prediction and the GW PSA in the 
Pn channel shown in FIG. [TT1 

As we see, the EFT in the strict heavy-baryon expansion works pretty well for observables 
(except for Pn), showing in the first two nontrivial orders good convergence to the data. However, 
there are also hints that the convergence is slow when one looks at how the parameters, particularly 
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fiA, change with order in TABLE [U It is also a bit surprising, but not necessarily significant, that 
r(<5) is found at NNLO to be much larger than Tew- 

The strict heavy-baryon expansion is in powers of Q/Mqcd and <5/Mqcd> and when Q ~ 
6 neither is a particularly small ratio. One can get a sense for the size of the corrections by 
considering the case where we retain some of the relativistic Q/m^ and 5/ tun corrections to all 
orders, the approach we call semi-resummation. We have already discussed in Sec. II V Dl how the 
slow convergence of the 5/m,M expansion affects the width. We now turn to its effects in comparison 
to data. 

At LO and NLO in the semi-resummation, there is no change in phase shifts compared to 
the strict heavy-baryon expansion, and therefore also no change in the position of the delta pole; 
only the relationship between width and parameters changes, and thus the parameters come out 
different. At NNLO, the P33 phase shifts, and with them the pole position, change because the 
dependence of the width on the energy is slightly modified. In FIG. [T2l we compare the NNLO 
results (Eqs. (|1Q3|) and (|104p ) of the partial resummation of relativistic corrections with those of 
the strict heavy-baryon expansion, exhibited before in FIGS. [101 and [TTJ The semi-resummation 
fit parameters, delta pole position, and scattering volumes are labeled "semi" in TABLES HI ILT1 
and IIIIl respectively. In this case we estimate the errors in 5/MeV, Jia, and x to be within ~ ±3, 
±0.20, and ±0.02, respectively, of the values in TABLE E At NNLO T(S) = 155 MeV. 

There is a slight improvement in the P33 phase shifts above the resonance, and the width at 
resonance is closer to the Breit-Wigner value, without destroying the good fit at lower energies. 
The fitted parameters vary less with order than in the strict expansion, but at NNLO parameters in 
two schemes agree within 5%, except for the small parameter x, where they differ by ~ 15%. The 
Pi 3 phases move in the right direction for most energies. In the other channels, the situation is less 
positive. The P31 phases get significantly worse. Finally, there is improvement in the Pn channel, 
but sizable higher-order contributions are still needed to get anywhere close to the empirical phase 
shifts. 

There are thus regions where each approach gets a better fit to data than the other. It is 
important to note that the semi-resummed curves are more or less within the error bands of their 
strict heavy-baryon counterparts, as one would expect from the fact that the differences between 
two methods are higher-order effects. Overall, we do not find that either the semi-resummation or 
the strict heavy-baryon works considerably better in describing data than the other. 
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FIG. 12: Comparison of NNLO P-wave phase shifts (in degrees) as functions of Wcm ( m MeV), the CM 
energy including the nucleon mass. The black solid (maroon dot-dashed) lines are the strict heavy-baryon 
(semi-resummed) results. The golden bands outline the estimated systematic errors of the semi-resummed 
expansion. The green dots are the results of the GW phase-shift analysis 



VII. SUMMARY AND CONCLUSION 



We have shown, using nN scattering as an example, how to generalize standard ChPT so that 
one can cope with the non-perturbative delta resonance within EFT. Our method is similar in 
spirit to that developed in Refs. 2CH22t|. but differs in detail. It can be thought of as a partial 



supporting argument for previous, successful unitarization results 



24 



25]. 



ChPT 

n 



ias been generalized to include an explicit field for the delta isobar. We built on earlier 



work [r], 3] based on the heavy-baryon formalism, which treats the delta as a nonrelativistic 
particle from the very beginning, rather than relying on a relativistic Lagrangian of the Rarita- 
Schwinger field. We worked out the 0(Q 2 /m 2 N ) relativistic corrections to the 7riVA vertex within 
this approach. 

EFTs are based on expansions in powers of the ratio of low to high scales, Mi /Mhi- The 
rationale to include the delta as an explicit degree of freedom is that an expansion makes sense 
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when M\ Q ~ m n ~ 5 and Mhi ~ Mqcd- Such an 



of phase space where the delta goes on shell 



expansion is straightforward away from regions 



13H15I]. but it requires change otherwise: certain 



contributions are enhanced when the external energy is dialed to the delta-nucleon mass difference. 
A way to deal with these enhancements had been explored in Refs. 2014221] . in which m y and 5 



were considered separate scales. Here we followed a similar approach to nuclear resonances [26] and 
constructed a power counting that incorporates the kinematic fine-tuning at the delta pole while 
keeping M\ Q ~ m n ~ 5. Our approach resulted in a global EFT description for both the threshold 
and resonance regions. 

Like other EFTs that deal with non-perturbative phenomena, ours captures the non- 
perturbative structure in LO. Subsequently, the power counting leads to a systematic, perturbative 
improvement beyond LO. We applied this power counting to low-energy ttN scattering, where we 
built the amplitudes up to NNLO. We have considered both a strict heavy-baryon expansion and 
a partial resummation of relativistic effects, which yielded comparable results. 

Up to NNLO, we could cast our P33 amplitude in terms of a Breit-Wigner form with energy- 
dependent width and background. Our EFT result is, however, not a mere reproduction of Eq. ([1]): 
with our approach, (1) the size of the background can be estimated beforehand; (2) one need not 
"cautiously" choose the proper energy domain for the resonance window; and (3) there is a link 
between the energy dependences in the width and background via the parameters 5, Ha, and x, 
which are then constrained by both the threshold and resonance data. 

We fitted our P-wave amplitudes to the phase shifts given by Ref. 28j |. With just three free 



pa_, we obtained a g „„d fit in the P33 channel. Contrary to an a, Hoc _tion in the 
resonance region [13], we had no difficulty preserving the unitarity condition around the resonance. 
We found parameters consistent with other determinations. With them, the P13 and P31 channels 
come out qualitatively correct. In the P\\ channel a sizable discrepancy exists between the EFT 
prediction, which is repulsive, and the empirical phase shifts, which become attractive within the 
region of study. 

In order to improve the description of the ttN data it is imperative to include other effects. One 
can push the calculation to next order (0(Q 2 )), when in channels other than P33 our framework 
should yield results similar to Ref. [3]. These lead to a better description of the phases, at least 
near threshold, and in particular were found to alleviate the P\\ discrepancy. The lack of attraction 
in this channel has been noticed long ago in the context of models, and it has been attributed to 
the Roper resonance's fairly large width (see, e.g., Ref. [3]). At 0(Q 2 ) Roper effects can appear 
through low-energy constants. Alternatively, Ref. 1?}] made an attempt of incorporating an explicit 
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Roper field in EFT but only investigated the region that is below the delta resonance. One can 
use the framework presented here to improve the description of the P\\ channel up to the energy 
of the Roper resonance. 

We have thus extended ChPT in ttN scattering to the delta region. The EFT approach pre- 
sented in this paper also provides the basis for a model-independent, unified description, from 
threshold to past the delta resonance without discontinuity, of reactions involving other probes 
and targets, including nuclei. These reactions, some of which have already been successfully stud- 
ied with the pioneering power counting of Ref. [^(J, present further tests of our power counting. A 
comprehensive confrontation of results from the two approaches should indicate which one provides 
a more efficient organization of EFT interactions. 



Acknowledgment s 

We thank Daniel Phillips, Vladimir Pascalutsa, James Friar, and (at a very early stage) Nadia 
Fettes for useful discussions. Comments on the manuscript by Daniel Phillips were greatly ap- 
preciated. We are grateful to the following institutions for hospitality during the long gestation 
of this work: Los Alamos National Laboratory (BwL), the Kernfysisch Versneller Instituut at Ri- 
jksuniversiteit Groningen (UvK), the National Institute for Nuclear Theory at the University of 
Washington (BwL, UvK), and the University of Arizona (BwL). This work was supported by the 
US DOE (BwL, UvK) and the Alfred P. Sloan Foundation (UvK). 



Appendix: Slow-velocity boosts 

A "heavy" particle is one that has three-momentum Q m, with m the mass of the particle. 
What we are looking for here is a Lorentz transformation rule for a heavy-particle field that is 
expressed in a perturbative fashion in powers of Q/m <C 1. For Q/m to remain small, boosts of 
the frame under consideration have to be in small velocities ~ Q/m. What we are studying here 
is in fact a special case of the heavy-baryon formalism, in which the four-velocity label v = (1, 0). 
Since v = (1,0) is sufficient for all the processes considered in this paper, we are not concerned 
with the invariance of the Lagrangian under a variation of v , namely, reparameterization invariance 



421 ] . More details will appear in Ref. [3^ 

A Poincare transformation takes a spacetime point x^ to 



c ,f * = A£x v + a* , (A.l) 
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with a four-vector representing the spacetime translation and the Lorentz-transformation 
matrix. The Poincare group has ten generators: the time translation H, spatial translations P, 
spatial rotations J, and boosts K. The commutation relations among these generators, or the 
Poincare algebra, can be readily found in the literature. We adopt the notation used in Ref. 331 ] . 



The most commonly used Poincare representations are a class of fields <&j(cc), with I a discrete 
index, that transform under the Poincare group as 

^ v (x) = M(A) in ^(x") , (A.2) 

where M(A) is a finite-dimension spacetime-independent matrix that furnishes a representation of 
the homogeneous Lorentz group and 

x" = A-^x-a) . (A.3) 

Therefore, for any (non-unitary) finite-dimension representation of the homogeneous Lorentz group, 
a corresponding Poincare representation can be constructed using Eq. (|A.2|) . Such representations 
include the well-known Klein-Gordon scalar field, Dirac field, four-vector field, etc. 

However, the Foldy-Wouthuysen (FW) representation [333 does not fall in the above category. 
A (classical) FW field with mass m and spin s, xi^^x), is an SO(3) (2s + l)-component spinor. 
The spin operators, s, are three (2s + 1) x (2s + 1) matrices, satisfying 

[si,Sj] = ie ijk s k . (A.4) 

The FW representation is furnished by identifying H, P, J, and K as follows, 

H = u , (A.5) 
P = p, (A.6) 
J = xxp + s, (A.7) 

- * I S X J) 

K = - (xu + ojx) — - tp, (A.8) 

with p = — iV and u = (m 2 + p 2 )z . One can check that the above definition does satisfy the 
Poincare algebra. In particular, an infinitesimal boost can be written as 

X '(t,x) = (l-ii-K)x(t,x) , (A.9) 

where the boost is parameterized by £, which is sometimes referred to as rapidity and is related to 
the relative velocity V between the original and boosted frame through 

£ = Vtanh^V . (A.10) 
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Because a; is a non-local operator, the FW representation is not as useful in relativistic situations 
as the Dirac field and the like. But the formal expansion of uj and K in powers of \p\/m enables 
one to derive a perturbative Lorentz transformation rule for the FW field, 

" = " i+ £-iS+-- i < A - u > 

K = mx — tp + - — (p 2 x + xp 2 ) sxp + -- - . (A. 12) 

4m v '2m 

We define the LO Lorentz transformation so as to reproduce the Galilean transformation. We take 
\p\ ~ Q, so that the kinetic energy E = oj — m scales as ~ Q 2 /m. Using the uncertainty principle, 
we expect that \x\ ~ l/Q and t~ 1/E ~ m/Q 2 . Also, we power count s as |s| ~ 1. Therefore, the 
LO boost generators are 

K(°) = mx - tp , (A. 13) 

and £ ~ Q/m. It is worth noting that does not depend on the spin operator, as indeed it 

should not. 

To remove the inert rest energy m, we define the heavy field 

^(t,x) = e imt X (t,x), (A.14) 

for which 

id t V(t, x) = {uj- m)*(t, x) . (A.15) 
An infinitesimal Galilean transformation of *(t,x) is found to be 

V\t,x) = (l-img-x + tg- V) *(t,f) = { 1 - imf -x + 0(Q/m)] ^} (*",£") , (A.16) 
where 

t" = (A- 1 x) Q = t + i-x , (A.17) 



X i 



{hr l x) 1 = Xi + tii . (A.18) 
For simplicity, we denote the boost transformation by 



l-im£-x * , (A.19) 



with the understanding that the left-hand side is evaluated at (t, x) while the right-hand side at 
{t",x"). 
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Going further down in the expansion in Eq. (|A.12|) . one finds the NLO boost generators, 



-> i 1 f 

= - — (p 2 x + xp 2 ) s x p , 

4m 2m 



(A.20) 



which depend on spin. The transformation rule for the nucleon field is worked out with s = <?/2, 

1 z* 



N 



l-im N £-x + - £-V+- £• K x v +°W K) 



N 



(A.21) 



2m at AmN 

When dealing with the delta, one needs to keep in mind that the mass-splitting 5 ~ Q in addition 
to the fact that the delta has spin 3/2, 



1 - im N £- x + — ^— C ■ V + i ■ (S& xv)+ 0(Q 3 /m 

2mjy 2m n V ' 



(A.22) 



Although 5 does not appear in the boost at this order, it does show up in the expansion of the 
kinetic energy. 

Our basic strategy to build an order-by-order Lorentz-invariant Lagrangian is the following: (i) 
enumerate all the rotationally invariant operators according to a certain scheme, e.g., the chiral 
index v defined in Eq. f)13[) : (ii) use Eqs. ()A.2ip and (jA.22j) to constrain the coefficients of these 
operators so as to preserve Lorentz invariance up to the order under consideration. Up to index 
v = 2 we find the terms shown in Eqs. (j!4|) . f)15[) . and f)16|) . One can explicitly check that the irNA 
operators, for example, are Lorentz invariant in the sense that Lorentz-violating effects only arise 
at v = 3 or higher. 
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